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(57) Abstract: The invention relates to a method of determining the complex impedance Z(f m ) of a non-steady electrochemical 
^ system, comprising the following steps consisting in: bringing the system to a selected voltage condition and applying a sinusoidal 
^> signal with frequency i m thereto; immediately thereafter, measuring successive values for voltage and current at regular time in- 
^ tervals ?T; calculating the discrete Fourier transforms for voltage (E(f)) and current (1(0), the voltage transform being calculated 
^ for the single frequency f m of the sinusoidal signal and the current transform being calculated for frequency f m and for two adje- 
ct cent frequencies f m .! and f m+1 on either side of frequency f m ; and calculating the impedance using the following formula: Z(f m ) 

0= E(f m yi*(f m ), wherein I* denotes a corrected current such that Re[I*(fJ] = Re[I(fJ] - |Re[Kf tn+1 )]+Re[I(f m .i)]}/2, Im[I*(f m )] = 
■ 



ImtKfJl - {ImlKUOl+ImrK^.Ol}^. 
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avani V expiration du dilai privu pour la modification des 
revindications, sera republiee si des modifications sont re- 
cues 



En ce qui conceme les codes a deux lettres et autres abrSvia- 
tions, se referer aux "Notes explicatives relatives aux codes et 
ubrevialions" figurunt au debut de chaque numero ordinaire de 
la Gazette du PCT. 



(57) Abrege : V invention concerne un precede de determination de l'impedance complexe Z(f m ) d'un systeme electrochimique non 
stationnaire, comprenant les etapes suivantes: niettre le systeme dans un etat de tension clioisi et lui appliquer un signal sinusoidal 
de frequence f m , mesurer immediatement apres des valeurs successives de la tension et du courant a des intervalles de temps regu- 
liers AT, calculer les transformees de Fourier discretes de la tension (E(f)) et du courant (1(f)), la trans formee de la tension etant 
calculee pour la seule frequence f m du signal sinusoidal et la transformee du courant 6tant calculee pour la frequence f m et pour deux 
frequences adjacentes f m _! et f m+ i de part et d'autre de la frequence f m , et calcucler l'impedance selon la formule suivante: Z^fJ = 
E(f m )/I*(f 1IL ) ou I* designe un courant corrige tel que: Re[I*(f m )] = Re[I(f ro )] - {Re[I(f ra+l )]+ReP(f IlM )]}/2, Im[I*(f m )] = Im[I(f m )] - 
llmrKUOHImrKf^,)])^. 
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PRO CEDE DE DETERMINATION DE L f IMPEDANCE D'UN SYSTEME 

ELECTROCHIMIQUE 

Pour la determination des proprietes et des qualites 
d'un sy steme electrochimique (tel qu'une pile, une batterie, un 
systeme de depot electrolytique , un systeme d» analyse d f un 
milieu) et la prevision de son fonctionnement futur, l'un des 
5 parametres dont la connaissance est essentielle est 1 1 impedance 

de ce systeme. 

De facon generale, un systeme electrochimique 
comprend, dans un milieu electrolytique, deux electrodes princi- 
pales - une electrode de travail et une contre- electrode. Une 

10 electrode de reference est disposee au voisinage de 1' electrode 
de travail et est utilisee en relation avec celle-ci pour 
effectuer diverses operations de mesure ou de regulation. Dans 
la suite du present texte les tensions et courants mentionnes 
correspondent a des mesures effectuees entre une electrode de 

15 travail et 1* electrode de reference correspondante . 

L 1 impedance d 1 un systeme electrochimique correspond au 
rapport tension sur courant. II est connu que pour mesurer rapi- 
dement la valeur de cette impedance en" fonction de la frequence, 
il est souhaitable de calculer le rapport des transformees de 

2 0 Fourier discretes (DFT) de la tension et du courant. Pour cela, 
on applique entre une electrode de reference et une electrode de 
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travail d'une cellule electrochimique un signal d' excitation de 
faible amplitude et on mesure a des intervalles de temps egaux 
AT, N valeurs successives e(n) et i(n) de la tension et du 
courant, avec 1 < n < N. L 1 expression generale des transf ormees 
5 de Fourier discretes E(f m ) et I (f m ) de la tension et du courant 
pour une sequence de N points est : 

N 

E(f m ) = AT £ etnJexp^TtjfnpAT) (1) 
n=l 

I(f m ) = AT £ i(n)exp(-27cjf IT pAT) (2) 

o€t j designe le nombre complexe dont le carre est -1, N le 

10 nombre de points de mesure, AT l'intervalle d 1 echantillonnage . 
Le calcul de transformed de Fourier discrete peut £tre effectue 
pour N/2 valeurs de frequence f m , avec 0 < m < (N/2) -1, N etant 
un nombre pair, Ces N/2 frequences discretes sont regulierement 
reparties entre 0 et 1/2AT (0. . .l/mAT. . . 1/2 AT) . 

15 L ' impedance complexe pour une frequence f m donnee est 

alors egale a : 

Z(f m ) = E(f m )/I(f m ) . (3) 
Les methodes de mesure d 1 impedance utilisees dans la 
pratique, different principalement par le type de signal d'exci- 

20 tation du systeme : sinusoidal, multi- sinusoidal, bruit blanc, 
etc . Une excitation sinusoidale est de loin le procede le plus 
utilise car il s'avere etre le plus precis. Le proc§de consiste 
a intposer une regulation sinusoidale a la cellule electro- 
chimique avec un signal de petite amplitude et a enregistrer la 

25 reponse de la cellule en courant et en tension. Le rapport des 
transf ormees de Fourier de la tension et du courant a la 
frequence du sinus donnera la valeur de 1 1 impedance a cette 
frequence. Le spectre des frequences est balaye en modifiant la 
frequence d 1 excitation. 

3 0 Dans la theorie des systemes electriques, l" expression 

de 1 1 impedance (equation (3)) est correcte a condition que le 
systeme analyse soit lineaire et stationnaire et que des pertur- 
bations ne soient pas apportees par des phenomenes exterieurs au 
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systeme. Dans le cas des systemes electrochimiqu.es, le respect 
de ces conditions impose des precautions particulieres. La 
fonction de transfert d'un systeme electrochimique est general e~ 
ment non lineaire mais on peut la considerer comme lineaire sur 
une petite portion. C'est pour cela qu'on utilise un signal 
d' excitation de faible amplitude. En consequence, le rapport 
signal /bruit diminue et on doit augmenter le temps de mesure, 
c'est-a-dire le nombre de points de mesure pour eliminer le 
bruit par une integration de la reponse. En revanche, 1" allonge- 
ment du temps de mesure pose des problemes quand a la stationna- 
rite du systeme electrochimique. Dans beaucoup de cas reels, la 
condition de stationnarite ne peut pas §tre remplie. Les causes 
sont multiples : relaxation du potentiel, relaxation du courant, 
relaxation des concentrations. Dans ces conditions, si le temps 
de mesure est suf f isamment long pour que 1 ' on ait un bon rapport 
signal/bruit, le systeme n'est pas stationnaire pendant la duree 
de mesure et 1 1 impedance calculee n f a plus grande signification, 
notamment aux basses frequences (par example inferieures au 
hertz) ou les problemes de non- stationnarite sont particuliere- 
ment aigus. 

Divers procedes ont ete proposes pour resoudre ce 
probleme de non stationnarite des systemes electrochimiques, 
mais aucun n'a donne de solution satisfaisante. > 



Ainsi, un objet de la presente invention est de 
pr6voir un nouveau procede de calcul de I 1 impedance d'un systeme 
electrochimique permettant de s'abstraire des erreurs liees a la 
non stationnarite du systeme. 

Pour atteindre cet objet, la presente invention prevoit 
un procede de determination de 1' impedance complexe Z(f m ) d'un 
systeme electrochimique non stationnaire, comprenant les etapes 
suivantes : 

mettre le systeme dans un etat de tension choisi et lui 
appliquer un signal sinusoidal de frequence f m/ 
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mesurer immediatement apres des valeurs successives de 
la tension et du courant a des intervalles de temps reguliers 
AT, 

calculer les transf ormees de Fourier discretes de la 
5 tension (E(f)) et du courant (1(f)), la transf ormee de la 
tension etant calculee pour la seule frequence f m du signal 
sinusoidal et la transf ormee du courant etant calculee pour la 
frequence f m et pour deux frequences adjacentes f m -i et f m +i de 
part et d 1 autre de la frequence f m , et 
10 calculer 1' impedance selon la formule suivante : 

Z(f m ) = E(f m )/I*(f m ) 
ou I* designe un courant corrige tel que : 

Re[I*(f m )] =Re[I{f m )] - {Re[I(f m+1 )]+Re[I(f m » 1 )] }/2 
Im[I*(f m )] = Im[I(f m )] - {Hn[I(f nH . 1 )]+Ita[I(f m - 1 )] J/2- 
15 Un mode de realisation particulier de la presente 
invention sera expose a titre non-limitatif en . relation avec la 
figure jointe qui represent e le spectre d l amplitude en fonction 
de la frequence d'un systeme electrochimique soumis a une 
excitation sinusoidale de faible amplitude et a un echelon de 
2 0 tension . 

Tout d'abord, on notera que, comme le potentiel 
applique est impose par l 1 instrument d' analyse, la non station- 
narite du systeme ne peut se manif ester que par des fluctuations 
de courant. 

2 5 La presente invention se base sur 1" analyse du compor- 

tement d'un systeme electrochimique soumis a un echelon de 
tension. La DFT de la reponse en courant traduit pour 
l ! essentiel des phenomenes de relaxation et done l ! effet de la 
non stationnarite du systeme. Comme le represente la courbe en 

3 0 point illes sur la figure unique, le spectre en amplitude I du 

module de la DFT du courant presente habituellement une forte 
reponse a basse frequence a des echelons de tension. 

Par contre, la reponse a une excitation sinusoidale de 
frequence f m d'un systeme stationnaire se traduit par une raie 
3 5 unique a la frequence f m . En pratique, cette reponse correspond 
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d'une part a la reaction du systeme a 1' excitation a la 
frequence f m , d 1 autre part a la contribution des effets de 
relaxation. 

L'avantage de 1 'utilisation de la transformee de 
5 Fourier est que le spectre monofrequence lie a l 1 excitation 
sinusoidale est superpose au spectre lie a. 1' Echelon de tension. 

Selon la pr£sente invention, une fois les considera- 
tions ci-dessus prises en compte, on propose de soustraire de la 
reponse en intensite a la frequence f m la contribution due aux 
10 relaxations, 6valuee a partir de 1' analyse de la DFT aux 
frequences adjacentes a f m . En effet, la transformee de Fourier 
pour des frequences f m -i et fm+1 ^ u systeme ne correspondra 
qu'aux phenomenes de non stationnarite et on considerera que la 
valeur de 1 ' amplitude liee aux non stationnarites pour une 
15 frequence f m est la moyenne des valeurs pour les deux frequences 
adjacentes. 

Ainsi, la presente invention propose : 

- de placer un systeme electrochimique dans un £tat de 
tension choisi et de lui appliquer une excitation sinusoidale 
2 0 monofrequence de faible amplitude, 

de mesurer, sans attendre la stabilisation du 
systeme, 1' amplitude du courant et de la tension a des interval - 
les de temps reguliers immediat ement apres 1 1 application de la 
tension tandis que l 1 excitation sinusoidale est appliquee, et 

2 5 -de calculer d'une part la transformee de Fourier 

discrete de la tension pour la frequence f m de 1' excitation et 
d' autre part la DFT du courant pour la valeur f m et pour deux 
valeurs de frequence f m -i et f m+ i adjacentes a la frequence f m , 
et 

3 0 -de calculer la valeur de l 1 impedance complexe a la 

frequence f m a partir de la valeur de la TFD de la tension pour 
la frequence f m et d'une valeur corrigee de la TFD du courant 
tenant compte des TFD calculees pour les frequences f m _i, f m et 

f m+l* 
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On calculera l 1 amplitude corrigee I* de la transf ormee 
de Fourier discrete du courant en valeur reelle et en valeur 
imaginaire par les deux equations suivantes : 

Re[I*(f m )] = Re[I(f m )] - {Re[I(f m+1 )]+Re[I(f m _ 1 )]}/2 
5 Im[I*(f m )] « Im[I(f m )] - {lm[l(f m+1 )]+lmtl(f m _i)] }/2 

La valeur de 1 ' impedance corrigee des effets de non 
stationnarite est alors : 

Z(f m ) = E(f m )/I*(f m ) 
avec I*(f m ) = Re[I*(f m )] + jlm[l*(f m )] 
10 Des le releve des points effectue pour la frequence 

f m , on pourra appliquer un signal sinusoidal a une nouvelle 
frequence et effectuer un nouveau releve, et ainsi de suite. 

Un avantage du present procede est qu'il rend possible 
une analyse correcte d 1 impedance sur des systemes electro- 
15 chimiques non stationnaires notamment aux tres basses 
frequences. En m§me temps, un gain de temps considerable est 
obtenu pour des systemes qui se stabilisent lentement car avec 
la correction selon 1' invention, il n'est plus necessaire 
d'attendre la stabilisation apres une mise sous tension du 
2 0 systeme pour demarrer une analyse d 1 impedance . 
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REVENDICATIONS 

1 . Procede de determination de 1 1 impedance complexe 
25 (f m ) d'un systeme electrochimique non stationnaire, caracterise 
en ce qu'il comprend les etapes suivantes : 

mettre le systeme dans un etat de tension choisi et lui 
appliquer un signal sinusoidal de frequence f m# 

mesurer ijnmediatement apres des valeurs successives de 
la tension et du courant a des intervalles de temps reguliers 
AT, 

calculer les transformers de Fourier discretes de la 
tension (E(f)) et du courant (1(f)), la transformee de la 
tension etant calculee pour la seule frequence f m du signal 
sinusoidal et la transformee du courant etant calculee pour la 
frequence f m et pour deux frequences adjacentes f m _i et f m+ i de 
part et d f autre de la frequence f m , et 

calculer 1 ■ impedance selon la formule suivante : 
Z(f m ) = E(f m )/I*(f m ) 
ou I* designe un courant corrige tel que : 

Ren*(f m )] - Re[I(f m )] - {RefKf^ll+RefKf^)]}^ 
Im[lMf m )] = Im[I(f m )3 - {lm[I (f m+1 ) ] +lm [I (%_!> ] }/2 . 

2. Procede selon la revendi cation 1, caracterise en ce 
qu'il est repete pour une succession de frequences d 1 excitation. 
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